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Abstract 

Background: Tuberculosis (TB) remains a major public health burden in many developing countries. China alone 
accounted for an estimated 12% of all incident TB cases worldwide in 2010. Several studies showed that the spatial 
distribution of TB was nonrandom and clustered. Thus, a spatial analysis was conducted with the aim to explore the 
spatial epidemiology of TB in Linyi City, which can provide guidance for formulating regional prevention and 
control strategies. 

Methods: The study was based on the reported cases of TB, between 2005 and 2010. 35,308 TB cases were 
geo-coded at the town level (n = 180). The spatial empirical Bayes smoothing, spatial autocorrelation and 
space-time scan statistic were used in this analysis. 

Results: Spatial distribution of TB in Linyi City from 2005 to 2010 was mapped at town level in the aspects of crude 
incidence, excess hazard and spatial smoothed incidence. The spatial distribution of TB was nonrandom and 
clustered with the significant Moran's / for each year. Local G* detected five significant spatial clusters for high 
incidence of TB. The space-time analysis identified one most likely cluster and nine secondary clusters for high 
incidence of TB. 

Conclusions: There is evidence for the existence of statistically significant TB clusters in Linyi City, China. The result 
of this study may assist health departments to develop a better preventive strategy and increase the public health 
intervention's effectiveness. 



Background 

Tuberculosis (TB) is an infectious disease caused by the 
bacillus Mycobacterium tuberculosis, which is still a 
leading cause of death in low-income and middle- 
income countries, and remains a major public health 
burden in many developing countries [1]. China (0.9 
million- 1.2 million) alone accounted for an estimated 
12% of all incident TB cases worldwide in 2010 [2]. Linyi 
City is one of the high-incidence areas in China, and 
there is an estimated 6,000 reported cases every year 
from 2005-2010. Thus, preventive strategies are ur- 
gently needed in the development of TB. 

In recent years, geographical information systems 
(GIS) and spatial analysis were frequently used to de- 
scribe the pattern of TB. In India, purely spatial and 
retrospective space-time analysis were used to find 
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significant hotspots of TB in three areas of the Almora 
district [3]. In Japan, space-time scan statistics identified 
TB clusters in Fukuoka [4]. In Madagascar, spatial clus- 
tering of TB was associated with socio-economic and pa- 
tient care factors in Antananarivo City [5]. With respect 
to spatiotemporal clustering, three high incidence space- 
time clusters were identified in Portugal between 2000 
and 2004 [6]. In South African, GIS and spatial analysis 
were used to study TB transmission patterns in a high- 
incidence area [7]. In an Urban West African, spatial 
scan statistic was used to assess purely spatial and 
space-time clusters of TB in Greater Banjul [8]. Another 
study in Brazil investigated spatial patterns of the inci- 
dence of TB and its relationship with socio-economic 
status [9]. In Beijing, GIS and spatial analysis were used 
to determine the role of migration in the transmission of 
TB [10]. To our best knowledge, there were few studies 
to explore the spatial epidemiology of TB in Linyi City, 
China. A better understanding of the spatial epidemi- 
ology of TB may help health departments to provide 
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guidance for formulating regional prevention and con- 
trol strategies [9,10]. 

The spatial analyses, such as spatial smoothing, spatial 
autocorrelation and cluster analysis are commonly used 
to characterize spatial epidemiology of diseases. Spatial 
smoothing was used to reduce random variation asso- 
ciated with small populations [9,11-14]. Spatial autocor- 
relation analysis was performed to detect significantly 
difference from a random spatial distribution [10,15-17]. 
Spatial cluster analysis was conducted to identify 
whether cases of disease were geographically clustered 
[4,8,14,18-25]. 

Therefore, we conducted GIS -based spatial analyses 
involving spatial smoothing, spatial autocorrelation ana- 
lysis and space-time scan statistic to characterize geo- 
graphic distribution pattern of TB in Linyi City of China 
during 2005-2010. 

Methods 

Data collection and management 

Relevant data for Linyi City (Figure 1), which is in the 
south of Shandong Province, was used in the study. 
Records on TB cases between 2005 and 2010 were 
obtained from the Linyi Institute for Tuberculosis Con- 
trol, and we were permitted to use the data. To conduct a 
GIS -based analysis on the spatial distribution of TB, the 
town-level polygon map at 1:100,000 scale was obtained, 
on which the town-level point layer containing informa- 
tion regarding latitudes and longitudes of central points of 
each county was created. Demographic information based 



on Linyi Statistical Yearbook was integrated in terms of 
the administrative code. All TB cases were geo-coded and 
matched to the town-level layers of polygon and point by 
administrative code using the software ArcGIS9.3 (ESRI 
Inc., Redlands, CA, USA). 

GIS mapping and smoothing 

To alleviate variations of incidence in small popula- 
tions and areas, annualized average incidences of TB 
per 100,000 at each town over the 6 year-period were 
calculated, and spatial empirical Bayes smoothing was 
implemented in SpaceStat software [26]. 

To assess the risk of TB in each town, an excess hazard 
map was produced. The excess hazard represents the 
ratio of the observed incidence at each town over the 
average incidence of all areas, the later was calculated by 
the number of cases over the total number of people at 
risk instead of the annualized incidence of a town [27] . 

Spatial autocorrelation analysis 

Spatial autocorrelation analyses were performed in Space- 
Stat. Global Moran s / statistics was used to discern spatial 
autocorrelation and detect the spatial distribution pattern 
of TB in Linyi City, China. Local G] was used to examine 
the local level of spatial autocorrelation and determine 
locations of clusters or hotspots. A calculated value of 
G\ >1.96 indicates that the town and its neighboring 
towns have a TB incidence rate statistically significantly 
higher than other towns. The number of permutation test 
was set to 999 and the significance level was set as 0.05. 




Figure 1 Location of the study area. 
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Space-time scan statistic 

The space-time scan statistic was performed using SaTS- 
can™ v9.Ll software [28]. The method is defined by a 
cylindrical window with a circular geographic base and 
with height corresponding to time [29]. The null hypoth- 
esis assumed that the relative risk (RR) of TB was the 
same within the window compared to outside. 

For this analysis, a Poisson based model was used, 
where the number of events in an area is Poisson dis- 
tributed according to a known underlying population at 
risk [30]. The geographic size of the window was limited 
to half the expected number of cases and that the time 
size was limited to half the total time period [29]. The 
test of significance of the identified clusters was based 
on comparing the likelihood ratio test statistics against a 
null distribution obtained from Monte Carlo Simulation 
[19]. The number of permutation was set to 999 and the 
significance level was set as 0.05. 



Results 

Descriptive analysis of TB in Linyi City 

There were 35,329 TB cases reported in Linyi City, from 
2005 to 2010. Of these, 35,308 (99.94%) had complete 
information including mapping of their place of resi- 
dence. Annualized average incidence at the town-level 
ranged from 27.21 to 100.58 per 100,000 (Figure 2). 

The excess hazard map showed distribution of the 
excess risk, which was defined as a ratio of the 
observed number over the expected number of cases. 
Towns in blue color had lower incidences than 
expected, as indicated by excess risk values less than 
1. In contrast, towns in red color had higher inci- 
dences than expected or excess risk values greater than 
1 (Figure 3). 

Spatial empirical Bayes smoothed map for annualized 
average incidence was created by correcting the variance 
in the variability of incidence (Figure 4). 
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Spatial autocorrelation analysis of TB in Linyi City 

The global spatial autocorrelation analyses for annual- 
ized incidence of TB in Linyi City from 2005 to 2010 
showed that the Moran s / was significant for each year 
(Table 1), implying that distribution of TB was spatially 
autocorrelated in Linyi City, China. 

Five significant spatial clusters (hotspots) of TB were 
identified using the Local G] for spatial autocorrelation 
(Table 2 and Figure 5). The hotspots persisted in Pingyi 
Sub- district, Gaozhuang Town, Xiawei Town, Mamuchi 
Town, Pingshang Town, Zhuanggang Town, Fangqian 
Town, Xiangdi Town, Zhubian Town, Shizilu Sub- 
district, Gaofengtou Town, Guichang Town, Lanshan 
Sub-district, Baishabu Town, Yitang Town, and Zao- 
goutou Town from 2005 to 2010. 

Space-time analysis of TB in Linyi City 

Space-time cluster analysis of TB in 2005-2010 in Linyi 
City showed that TB was not distributed randomly in 
space-time. The most likely statistically significant cluster 
for high incidence of TB was found to exist at Pingyi Sub- 
district, for the year 2008-2010 (RR = 2.22, p < 0.001), 



with 520 observed cases and 235.93 expected cases. 
Nine statistically significant secondary clusters were 
also detected for high incidence of TB. The results 
are listed in Table 3, and depicted on the map in 
Figure 6. 

Discussion 

In our study, exploratory spatial data analysis and spatial 
cluster analysis of TB were conducted at town level in 
Linyi City, China. We mapped TB from different aspects 
such as crude incidence, excess risk, and spatial empir- 
ical Bayes smoothed incidence, investigated the spatial 
pattern and highlighted geographic areas with significant 
high incidence of TB in Linyi. The study showed that 
the spatial distribution of TB in Linyi City was nonran- 
dom and clustered with the significant Morans / for 
each year. Local G] detected five significant spatial clus- 
ters for high incidence of TB when only space distribu- 
tion was considered. However, one most likely cluster 
and nine secondary clusters for high incidence of TB 
were identified when both space and time were consid- 
ered in the space-time analysis. When compared the 




0 15 30 45 60. 



Kilometers 

Figure 3 Excess hazard map of TB in Linyi City, China from 2005 to 2010. 
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Figure 4 Spatially smoothed map of TB in Linyi City, China from 2005 to 2010. 



clusters of the Local G\ with those of the space-time Table 2 Spatial clusters (hotspots) of the TB identified by 
scan statistic, both methods detected similar and signifi- local G f * in Linyi City, China from 2005 to 2010 



cant high-risk clustering. Consistent results using these 
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The result of the present study provided useful infor- 
mation on the prevailing epidemiological situation of TB 
in Linyi City. The novel knowledge about the presence 
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Table 1 Global spatial autocorrelation analyses for 
annualized incidence of TB in Linyi City from 2005 to 
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of clusters of TB in Linyi can help the Linyi Institute for 
Tuberculosis Control to intensify their remedial mea- 
sures in the identified areas of high tuberculosis preva- 
lence and chalk out future strategies for more effective 
TB control. Strategies may include compulsory BCG 
immunization of the children, educating the public 



about the dangers of TB as a re-emergent epidemic, and 
monitoring TB carefully. 

While our study has demonstrated the usefulness of 
GIS and spatial analysis, it still has several limitations. 
First, our data relies on official surveillance and we can- 
not exclude the possibility that some towns may 



Table 3 SaTScan statistics for space-time clusters with significant higher incidence in Linyi City, China from 2005 to 
2010 
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Figure 6 Space-time distribution of the detected clusters of TB cases with significant higher incidences in Linyi City, China from 2005 
to 2010. 



underreport the number of cases forvarious reasons. 
Cases might be missed by routine notification systems 
because people with TB do not seek care, seek care but 
remain undiagnosed, or are diagnosed by public and pri- 
vate providers that do not report cases to local or na- 
tional authorities. Second, we analyzed a relatively short 
period of time (i.e. 6 years, from 2005 to 2010). Further 
studies are needed to evaluate the spatial and temporal 
changes in the pattern of TB using data from a longer 
study period. Third, we did not assess possible risk fac- 
tors that could be associated with clustering. It is not a 
survey-based study but official surveillance, as socio- 
economic and environmental factors are not taken into 
account. 

The present study only analyzed the statistically 
significant clusters of TB. Future researches are war- 
ranted to focus on the effect of various socio-economic 
and environmental factors on the high incidence of TB 
in the clustering areas. Disease prevalence is frequently 
associated with many aspects of socio-economic status, 



such as overcrowding [7], unemployment [7,31], low 
educational level [32,33], number of shebeens [7] and 
poor housing quality [34]. Moreover, spatial clustering 
of TB was also associated with the migrant popula- 
tion [10], patient care factors [5] and environmental 
factors [7]. After detecting the statistically significant 
clusters of TB in the region, a survey-based study is 
intended to identify the role of these factors in the 
spread of TB. 

Conclusion 

This study showed the presence of clusters of TB in 
Linyi City, China and demonstrated that using the exist- 
ing health data, the spatial statistic and GIS may provide 
public health officials with necessary feedback about the 
prevalence of statistically significant clusters of TB in 
the region, and thus enable them to carry out more ef- 
fective strategies to control TB. Since the efficacy of TB 
control measures in specific areas could be assessed by a 
longitudinal change in TB prevalence, the space-time 
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scan statistic also can contribute to health program 
evaluation. More detailed individual level investigations 
are needed in the identified clusters to evaluate the most 
important determinants of disease distribution. 
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